Nonlinear optics of semiconductors under an intense terahertz field 
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A theory for nonlinear optics of semiconductors in the presence of an intense terahertz electric field 
is constructed based on the double-line Feynman diagrams, in which the nonperturbative effect of 
the intense terahertz field is fully taken into account through using the Floquet states as propagating 
lines in the Feynman diagrams. 
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I. INTRODUCTION 

Since early 1990's, thanks mainly to the emergence 
of free electron lasers operating in the terahertz (THz) 
waveband^ the interaction between semiconductors and 
a strong THz field has been brought under intensive 
investigations. Nonlinear transporlsSi&iiii and linear 
opticsS*2*2*2iiSiiiiiS are the two main themes of these in- 
vestigations (here we use the terminology of "linear op- 
tics" or "nonlinear optics" in the sense that the intense 
THz field is treated as a part of the system but not an 
optical excitation, otherwise, if the THz field is viewed as 
an external optical field, even the so-called "linear optics" 
here would be highly nonlinear). To thoroughly under- 
stand the physics in THz-field-driven semiconductors, as 
well as to develop novel devices based on these systems, 
nonlinear optical spectroscopies are a powerful and some- 
times necessary method due to their accessibility both 
in ultrafast time-resolution and in multi-frequency mix- 
ing. For example, the four-wave mixing spectroscopy 
has been adopted to study the effect of the strong THz 
field induced by Bloch oscillation in biased semiconductor 
superlatticesyi^ and a theory based on Floquet states^ 
of time-periodic systems has been developed to consider 
the non-perturbative effects of the THz dipole field.— Re- 
cently, the difference-frequency processes were proposed 
to generate THz emission, and the estimated strength of 
THz field could be of the order of kV/cmjiSiii which is 
so large that the feedback effect of the THz field on the 
nonlinear difference- frequency process may be important. 
So, as a common theoretical basis, a nonlinear optical 
theory of semiconductors in the presence of an intense 
THz field is desired. 

To construct such a theory, it is essential to include 
the non-perturbative driving of the THz field. To this 
end, a good starting point is the eigen states of the THz- 
driven systems, the Floquet states^ which have the non- 



perturbative effect of the THz field fully included. In 
fact, a compact theory for linear optics of THz-driven 
semiconductors has been formulated in the Floquet-state 
basis£ In next section, the general formalism for nonlin- 
ear optics of semiconductors under a strong THz field will 
be constructed with the double-line Feynman diagrams 
frequently used in textbooksiiS'i^ In Section IIIII some 
examples will be given to illustrate how to calculate the 
nonlinear optical susceptibility from the Feynman dia- 
grams. And the conclusions are given in the last section. 



II. GENERAL THEORY 

The system to be considered is a semiconductor irradi- 
ated by an intense cw THz laser. The Hamiltonian of this 
system under excitation of additional weak lasers can be 
expressed as 



IT = JT (t) - X> • F„(f), 



(1) 



where Ho(t) is the unperturbated Hamiltonian of the 
semiconductor with the THz-field-driving included, fi is 
the dipole operator, and 
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is the pertubative optical field. 

With the density matrix of the system denoted by p(t), 
the optical polarization is P(t) = Tr [p(t)fi\. As the THz 
field, with photon energy much smaller than the band 
gap, induces no inter-band excitation, the system is as- 
sumed in the semiconductor ground state before optical 
excitation, i.e. p(—oc) = |0)(0|. Thus the jth component 
of the Nth. order [^W] nonlinear optical response to the 
optical fields isiLiS 
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where the summation is over all permutations of the in- 
dices as indicated, and 

U{t,t') = Te-*St> Ho[tl)dtl = UHt',t) 

is the unperturbed propagator of the system. The sys- 
tem in the presence of an THz field is time-periodic, i.e. 
H (t) = H {t + T), where T = 2tt/uj with uj denoting 
the angular frequency of the THz field. The eigen states 
of the time-periodic Hamiltonian is the Floquet states 
{|<Z,t)}, which are time-periodic and satisfy the secular 
equation 

[H (t) -ihd t ] \q,t) = E q \q,t) = E q \q,t + T), (4) 

where E q , in analogy of quasi- momentum of Bloch states, 
is termed quasi-energy. Obviously, the sidebands of the 
Floquet states \q,m,t) = cxp(imujt)\q,t) are also eigen 
states of Eq. Q with quasi-energy E qiJn = E q + mfiuj. 
Now the propagator can be expanded into the Floquet 
states as 



U(t,t') = \q,t)e-K E ^- t 'Hq,t , \, 



(5) 



(hereafter all superscripts and subscripts appearing only 
on the righthand side of an equation are assumed dumb 
indices to be summed over). The dipole matrix element 
between the Floquet states is 



»q;q'(t) = (q,t\fl\q' ,t) 



imuit 



(6) 



\q,n><q,n+m\ 

A 



\j,m> 



£2, 



(a) 



<c,m I 

a 

c 

<a,m I 



10x01 



((1) IX^WjXXpWj+ml 
D. 



A A 
t 



<0,n\ 



<X 2 ,n 2 \ 



ft, 



10x01 



\X,nxX,n+m\ 



£1, 



(b) 



A A 
t 

h 

10x01 



Q. 



10x0, ml a 



A A 
t 

t 



<X,n\ 
(C) 10x01 fi i 



(e) \B,n><B,n+m\ 



A A 
t 



\X v n t > 



ft, 



<X y n 3 \ 



10x01 



where 



(q, m, t\fi\q', t)dt = fi q . q > _ T 



is the time-average of the dipole matrix element. Thus 
we have 

U(t',t)fr = \q',m\t')e-* E «'^'<- t '~V(i ql , m ,; q (q,tl(7) 



ftU{t,t') = \q,t)n q ; 



q' .rn' 



(q\ m', f | . (8) 



With Eq. (0) and © used respectively on the left 
and right sides of |0}(0| in Eq. ©, the nonlinear optical 
response can be derived. The double-line Feynman dia- 
grams can be used to assist such a derivation as well as in 
static cases. For this particular time-periodic system, the 
rules of the Feynman diagrams for the optical response 
P^ N) (t) are [see Fig. [2(a)]: 

1. The evolution of the diagram starts from the 
ground state |0)(0| at t — — oo and ends at t with 
the final Floquet-state density matrix \q,n)(q,n + 
ml. 



FIG. 1: (a) The general form of double-line Feynman dia- 
grams for nonlinear optics of semicondcutors in the presence 
of an intense THz electric field, (b), (c) The double- line Feyn- 
man diagrams for linear optics of THz-field-driven semicond- 
cutors. (d),(e) Two of the double-line Feynman diagrams for 
X*- 3 - 1 four-wave mixing in THz-field-driven semiconductors. 



The N interaction vertices at t p < t (p = 
1,2, ... , N) consist of a photon line (the dotted ar- 
row) with frequency f2 p , the Floquet state before 
interaction (|a,m a ) on the left branch or (a, m a \ 
on the right branch), and the state after interac- 
tion (|c, m c ) on the left branch or (c, m c | on the 
right branch). The photon lines pointing to inside 
and outside represent photon absorption and emis- 
sion processes, respectively. The pth vertex con- 
tributes a dipole matrix element j-H c ,m c ;a,m a (or 
"fl J a,m a ;c,mJ if it's on the left (or right) branch, 
and a factor from the optical field F p exp(— ifl p t p + 
iK p • R) for photon absorption or F* exp(ifl p t p — 
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iK p ■ R) for photon emission. 

3. The photon emission vertex at final time t con- 
tributes the factor (fij) c ,m a -q,n+me~ lmult where 
(c, m c \ is the state before photon emission. 

4. The double-line between two neighbor vertices at t\ 
and ti [t\ < t%) represents the unperturbed propa- 
gator of the density matrix. If the state associated 
with the double- line is \b, rrib}(c, m c \, the factor due 
to the propagator is 

9(t 2 - t 1 ) e -*( B »"»»' -™»- <r »-»»> -"»»)( t!, -* 1 J, 

where E^ rrlb . c , mc = E h<mh - E c , mc is the tran- 
sition energy between the Floquet states, and 
^b,m b -c,m /fi is the relaxation rate of the density 
matrix due to the interaction with environment. 



The rules for constructing susceptibility in frequency do- 
main from the Feynman diagrams can be easily derived 
from the rules above: 

1. The interaction vertex, consitituted by the pho- 
ton line with frequency Q p , the initial state 
\a,m a ) (or (a,m a \), and the final state |c, m c ) (or 
(c,m c \), contributes the factor j-(fJ-j p ) c .m c ;a,m a (or 
_ l(A t ip)a,m a ; C ,m c ) if the vertex is on the left (or 
right) branch . 

2. The factor associated with the double-line state 
\b, m/,)(c, m c \ between t\ and £2 (^i < ^2) is 

t n <t 2 

3. The factor associated with the final photon emis- 
sion vertex at t is 

N 
n=l 

where (c, m c | is the state before photon emission. 

4. All factors are multiplied together and all inter- 
mediate Floquet states and their sidebands are 
summed over. Then all possible diagrams are 
summed. For a certain optical configuration, the 
frequency Q p is positive or negative depending on 
whether the corresponding vertex is an absorption 
or emission process. 



5. For each diagram, all factors described above 
should be multiplied together and all Floquet states 
and their sidebands are summed over. Then all dif- 
ferent diagrams, which are determined both by how 
the N vertices are grouped into two branches and 
by the time-ordering, should be summed. For a cer- 
tain optical configuration (which determines each 
vertex is a photon absorption or emission process), 
there are totally 2 N N\ different diagrams, among 
which there are 2 N elementary diagrams and the 
others can be derived from them by permutation of 
the N vertices. And still, there are 2 N different op- 
tical configurations. In a certain experiment, only 
one optical configuration needs to be considered. 

The nonlinear optical susceptibility in frequency do- 
main is defined as 



(9) 

I 

From the rules described above, we can see that the 
nonlinear optical response of THz-field-driven semicon- 
ductors, or generally speaking, time-periodic systems, 
takes the form very similar to the textbook formalism for 
static systemsiiSiiS The difference lies on three aspects: 
First, the dipole matrix element here is the time-average 
of that between Floquet states. Secondly, the THz field 
induces new resonances at THz-photon-assisted transi- 
tions, as the sidebands of the Floquet states act as in- 
termediate states in the nonlinear optical process . And 
thirdly, there is an extra dynamic phase factor e - imuJt at 
the final photon-emission vertex, which makes the energy 
of the emitted photon differ from the total input energy 
by an integer multiple of the THz-photon energy, corre- 
sponding to the physical process of THz-photon-sideband 
generationi&i2i In the time-periodic systems, the energy- 
conserving condition is relaxed to the quasi-energy con- 
servation. We reiterate that the nonperturbative effect of 
the THz field has been fully included through the renor- 
malization of dipole matrix element and transition en- 
ergy. 



III. EXAMPLES 

As an illustrative example, the Feynman diagrams for 
linear optics are plotted in Fig. (b) and (c). From the 
rules for the Feynman diagrams, the linear susceptibility 
can be formulated as 



, jjv)e 
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^1 Eo-.X.n ~ «ro ; X,i 



(10) 



where |X, n) denotes the THz-photon sidebands of the 
Floquet-state excitons and \0,n) denotes the sidebands 
of the ground state . This result is identical to that in 
Ref. |8j derived with non-equilibrium Green's function 
technique. 

Now we consider another example, the four-wave 
mixing, in which three input beams propagate in the di- 
rections Ki, K2, and K3, respectively, and the signal is 



detected in the direction Ki + K2 — K3. Corresponding 
to this optical configuration, there are 48 different Feyn- 
man diagrams, among which only 16 are under resonant 
excitation condition. We calculate two typical diagrams 
as examples of resonant excitation of excitons and bi- 
excitons (exciton molecules constituted by two excitons) . 
Fig. ^ (d) is a Feynman diagram for resonant excitation 
of excitons, which contributes to the susceptibility as 



27n5(» -Q 1 -n 2 +n 3 - mu) (Mj3)o;X 2 ,» 2 Mx 2 ,n 2 ;0,n (MjQo.njXx.m+m ^h) X u nyfi 

(Ex uni ;0,n — iTxi,m;0,n — f^l — ^2 + ^3) {Exr ,ru ;X 2 ,n 2 ~ X l .n l -X 2 ,n 2 — ^1 + ^3) (E 0; x 2 ,n 2 ~ i^Q;X 2 ,n 2 + ^3) ' 

and Fig. ^ ( e ) is a diagram for resonant excitation of bi-excitons, which contributes to the susceptibility as 

-27nS(n - r>i - n 2 + n 3 - mop (n j3 ) . X3in3 (^)x 3 ,n 3] B,n+ m i^h^Si^i ^x 1 , nii0 

(EB,n;X 3 ,n 3 ~ i^B,n;X 3 ,n 3 ~ ^1 — ^2 + ^3) (Exx ,m;X 3 ,n 3 ~ fi?Xi ,ni \X 3 ,n 3 — Ol + ^3) (E Xl . ni ;0 ~ ^Xi,ni;0 ~ ^l) ' 



where \B,n) denotes the sidebands of Floquet-state bi- 
excitons. From the two terms above, we can easily 
identify the resonances associated with the THz-photon- 
assisted Floquet-state exciton and bi-exciton transitions. 
The sideband generation, as indicated by the 5-function 
with argument containing an integer multiple of THz- 
photon energy, accounts for the four- wave mixing signal 
out of the excitation spectrum observed in the numerical 
calculations in Ref. 15. 



IV. SUMMARY 

In summary, based on the double-line Feynman dia- 
grams similar to the static case in textbooks, we have 



I 

construct a general theory for nonlinear optics of semi- 
conductors under an intense THz field. The basis of the 
Feynman diagram is the eigen states of the time-periodic 
systems, i.e., the Floquet states, so the non-perturbative 
effect of the THz field has been fully taken into account. 
Many phenomena, including the THz-photon- assisted ex- 
citon or bi-exciton resonances and the THz-photon side- 
band generation, are naturally accounted for in this the- 
ory. 
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